% 参数设置
k = 5; % 弹簧劲度系数
m = 20; % 两端质点的质量
mo = 20; % 中间质点的质量

% 初始条件
x0 = [0; 0; 0]; % 初始位移
v0 = [0; 0; 0]; % 初始速度

% 使用拉普拉斯变换法求解系统的运动微分方程，得到解析解
% 这部分需要根据具体的系统方程来实现
% 假设我们已经得到了解析解 x1(t), x2(t), x3(t)

% 时间参数
t = linspace(0, 10, 1000); % 时间跨度

% 计算位移
x1 = ...; % 根据解析解计算x1随时间的变化情况
x2 = ...; % 根据解析解计算x2随时间的变化情况
x3 = ...; % 根据解析解计算x3随时间的变化情况

% 创建动画
fig = figure;
ax = axes('Parent', fig);
line([0, x1(1)], [0, 0], 'Color', 'b');
line([x1(1), x2(1)], [0, 0], 'Color', 'r');
line([x2(1), x3(1)], [0, 0], 'Color', 'g');
xlim(ax, [-1, 4]);
ylim(ax, [-1, 1]);

for i = 1:length(t)
    set(gca